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Abstract. We show that the Cosmic Microwave Background (CMB) polarization data 
gathered by the BOOMERanG 2003 flight and WMAP provide an opportunity to investigate 
in-vacuo birefringence, of a type expected in some quantum pictures of space-time, with a 
sensitivity that extends even beyond the desired Planck-scale energy. In order to render this 
constraint more transparent we rely on a well studied phenomenological model of quantum- 
gravity-induced birefringence, in which one easily establishes that effects introduced at the 
Planck scale would amount to values of a dimensionless parameter, denoted by ^, with respect 
to the Planck energy which are roughly of order 1. By combining BOOMERanG and WMAP 
data we estimate ^ ~ — 0.110 ± 0.076 at the 68% c.l. Moreover, we forecast on the sensitivity 
to ^ achievable by future CMB polarization experiments (PLANCK, Spider, EPIC), which, in 
the absence of systematics, will be at the 1-a confidence of 8.5 x 10~^ (PLANCK), 6.1 x 10"'^ 
(Spider), and 1.0 x 10~^ (EPIC) respectively. The cosmic variance-limited sensitivity from 
CMB is 6.1 X 10^6. 
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1. Introduction 

The challenge of finding a theory for "quantum gravity" to reconcile quantum mechanics 
with the General Relativity description of gravitational phenomena has been confronting 
the theoretical physics community for more than 70 years [Ij. A reason why after many 
decades we understand little on the "quantum-gravity problem" originates from the difficulties 
encountered in experimentally accessing the realm of quantum corrections to gravity, expected 
to be originating at the ultrahigh energy scale given by the Planck energy, Ep ~ lO^'^eV (or 
equivalently the ultrashort Planck length, Lp = 1/ Ep 10~^^m). 

Over the last decade the search for experimental manifestations of Planck-scale effects 
has been reenergized by the realization that certain types of observations in astrophysics do 
provide indirect access to certain Planck-scale effects [21 El IH El [6], [7]. Several authors have 
already argued that cosmological observations may soon also help in the efforts of searching 
for hints on the realm of quantum gravity. These expectations originate from the fact that 
many cosmological observations reflect the properties of the Universe at very early times, when 
the typical energies of particles were signiflcantly closer to the Planck scale than the energies 
presently reached in our most advanced particle accelerators. Moreover, the particles studied 
in cosmology have typically travelled ultra- long ( "cosmological" ) distances, and therefore even 
when they are particles of relatively low energies they could be affected by a large accumulation 
of the effects of the "space-time quantization" , which is one of the most common expectations 
emerging from quantum-gravity research^. 

Over the last few years some authors (see, e.g., Refs. [HlEllinilll]) have indeed started to use 
ideas coming from quantum-gravity research in cosmology, but the possible connections between 
cosmological measurements and quantum-gravity research have yet to be fully appreciated in 
the literature. This is because in the relevant cosmological studies, while one indeed easily 
sees that the nature of a certain effect is of interest for constraining models of quantum 
gravity, it is often not easy to establish whether the conjectured magnitude of the change could 
plausibly be connected with corrections introduced genuinely at the Planck scale. Our main 
objective here is to establish that in the case of studies of quantum-gravity-induced in-vacuo 
birefringence, proposed in several quantum-gravity studies (see, e.g., Refs. [3l[12l[l3]) with the 
CMB polarization data we may indeed reach sufficient sensitivity to test effects introduced at 
the Planck energy scale. 

Several studies have already investigated (see, e.g., Refs. [HI UHl [lEl [171 UHl [I9]) the 
relevance of CMB polarization data for testing birefringence, and some studies have motivated 
the analysis also from a broadly- intended quantum-gravity perspective [TTl [20], but there is 
still no dedicated study attempting to establish whether the sensitivity levels provided by 
CMB observations could provide access to birefringence introduced at the Planck energy scale. 
Here, we make this connection clear and provide a constraint using existing BOOMERanG and 
WMAP polarization data. Moreover, we also make predictions for the expectations with future 
CMB polarization measurements both with space and sub-orbital experiments. 

In order to make our case for Planck-scale sensitivity more transparent we model 

§ Most approaches to the quantum-gravity problem predict some form of space-time quantization, although 
different approaches often lead to profoundly different predictions on what space-time quantization should entail. 
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birefringence in the context of a much-studied formahzation of Planck-scale modifications of 
electrodynamics, a model proposed by Myers and Pospelov |[T2] , with Lagrangian density of the 
form 

^ = -If.^F^" + ^f'PoAFui . (1) 

This model was proposed as a field theory formalization of the intuition on Planck-scale 
modifications of the energy-momentum dispersion relation that had been developed in a series 
of preceding studies starting from those in Refs. [21 [3l [21]. The Lagrangian has been much 
studied from the quantum-gravity perspective, including several analyses on its observable 
features (see, e.g. Refs. [3 [131 [221 [23l [211 [25]). 

We feel that the fact that our suggestion for testing Planck energy scale corrections in 
CMB polarization data is based on such a well-studied model should render our findings more 
significant. Moreover, in the case of the Myers- Pospelov field theory it is relatively simple 
to establish which range of values of the single parameter of the model, ^, would amount to 
introducing the effects at the Planck energy scale: the correction term in ([T|) is of dimension 
5, so that the parameter ^ must be dimensionless, and several studies [3, [121 ESI ESI [2ll [25] 
have shown that quantization of space time at the Planck scale would lead to correction terms 
with ^ roughly of order 1 (|^| ~ 1) in equation ([T]). This estimate of course must be handled 
prudently, since the arguments estimating that the scale of quantum-gravity effects should be 
given by the Planck scale allow a certain range of uncertainty: for example, most of these 
arguments rely crucially on the rather optimistic assumption of having no running of the 
gravitational coupling constant and no role for gravitational effects in the running of other 
coupling constants, but plausible estimates of how the renormalization-group flow might be 
affected by gravitational effects show [26] that the quantum-gravity scale may well differ 
by 2 or 3 orders of magnitude from the Planck scale. It is therefore rather significant that the 
analysis we report here establishes robustly that available CMB polarization data are sensitive 
to |,^| ~ 1 and that planned more refined measurements of the CMB polarization will provide 
access to values of |^| even a few orders of magnitude smaller than 1. 

Interestingly, through a combined analysis of BOOMERanG and WMAP data as we 
describe in detail later we find the following estimate: ^ ~ —0.110 ±0.076 at the 68% c.l.. This 
constraint very explicitly shows that indeed the relevant CMB polarization data have sensitivity 
at Planck-scale energies. And it is also intriguing that the analysis leads us to an estimate of 
^ which is nonzero at roughly the 1.5(7 level. We shall stress that most previous analyses of 
birefringence in relation to BOOMERanG and WMAP data assume that the polarization angle 
rotation, usually denoted by a, is energy independent, a = ao, whereas within the Myers- 
Pospelov framework one finds a characteristic energy dependence of a: a = ^E'^T, where E 
is the energy of the radiation and T is the time of propagatioijjj]. Using the same data that 
lead to our result ^ = — 0.110±0.076 one would instead find ao = (— 2.4±1.9)° assuming an 
energy independent rotation. One could therefore (however weakly) argue that the evidence in 
favor of birefringence is somewhat more robust for the case of the energy-dependent effect we 
have considered (a 1.5cr effect versus a 1.3o" effect). 

II This dependence is inverse to the energy dependence expected in polarization rotation associated with Faraday 
rotation in the presence of a magnetic field with the rotational angle proportional to E^^. 
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We also comment on previous attempts to constrain ^ using astrophysical data, such as the 
ones reported in Refs. [3 [221 ES]- These analyses provide bounds on that naively appear to 
be tighter than the ones we obtained here. However, we find that the frame dependence of the 
analysis of such symmetry-breaking effects renders those astrophysical analyses inapplicable to 
the context of CMB studies, performed of course in the reference frame naturally identified by 
the CMB. A crucial point for this observation is the realization that the space-isotropic form 
of the Myers-Pospelov Lagrangian density can at best be assumed in a very restricted class of 
frames of reference: even when the Lagrangian density takes space-isotropic form in one class of 
frames it will not be in general space-isotropic in other frames. And our reasoning leads us also 
to conclude that in the hypothetical scenario of future more sensitive CMB polarization data 
providing more robust evidence of energy-dependent birefringence (say at a 3a level), it would 
be inappropriate to rush to the conclusion that the space-isotropic Myers-Pospelov framework 
is finding confirmation. In such a scenario one should then necessarily perform searches of 
possible evidence of anisotropy. 

The paper is organized as follows: In the next Section, we describe the Myers-Pospelov 
framework, and its implications for birefringence. Then in Section 3 we discuss some aspects of 
the phenomenology of the Myers-Pospelov framework that are particularly relevant for studies 
of the CMB. Section 4 contains our key result obtained from CMB polarization observations. 
In Section 5 we compare our analysis to previous, partly related, analyses, and in particular 
we discuss the issue of frame-dependent spatial isotropy. We conclude, in section 6, with a 
brief summary of our results and some observations on the outlook of the studies we are here 
proposing. We shall work throughout adopting "natural units" {h = c = 1), and conventions 
rj^'^ = diag(l, —1, —1, —1), for the metric, and £123 = 1, e^'^^ = —1, for the Levi-Civita symbol. 

2. Effective field theory for Planck-scale-modified electrodynamics 

Motivated by previous studies [21 [3l EI] which had discussed some mechanisms for Planck- 
scale modifications of the laws of propagation of microscopic particles, in Ref. [12] Myers and 
Pospelov proposed to describe such effects within the framework of low-energy effective field 
theory, and observed that, assuming essentially only that the effects are linear in the (inverse 
of the) Planck scale and are characterized mainly by an external four-vector, one arrives at a 
single possible correction term for the Lagrangian density of electrodynamics: 

C = --F^.F^^'' + ^n''F^sn''d^{npe^''^^F,^) (2) 

where the four- vector parameterizes the effect. 

The (dimensionless) components of of course take different values in different reference 
frames, transforming indeed as the components of a four-vector. The quantum-gravity intuition 
that motivates the study of this model also inspires the choice of the factor and leads one to 
the expectation that the components of the Myers-Pospelov four-vector should (when nonzero) 
take values somewhere in the range 10"'^ < Ua < 10^. Values of order 1 would essentially 
correspond to introducing the effects exactly at the Planck scale. Values as high as 10^ would 
still be plausible, especially in hght of some renormalization-group arguments suggesting [26] 
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that the characteristic scale of quantum-gravity effects might actually be closer to the particle- 
physics "grand-unification scale" than to the Planck scale. There is no robust argument that 
would suggest that the characteristic scale of quantum-gravity effects might instead be much 
higher than the Planck scale, but most authors prudently consider values of parameters such as 
the n^'s still meaningful even down to 10~^, because of the desire to be rather cautious before 
excluding the possibility of a quantum-gravity interpretation. 

These arguments provide motivation for developing a phenomenology of the Myers- 
Pospelov field theory exploring a four-dimensional parameter, nc, keeping in focus the most 
interesting range of values, 10^^ < < 10^, and contemplating the characteristic frame 
dependence of the parameters tIq,. There is already a rather sizeable literature on this 
phenomenology (see, e.g., Refs. [211 E5] and references therein) but still fully focused on 
what turns out to be the simplest possibility for the Myers-Pospelov framework, which is the one 
of assuming to be in a reference frame where ria only has a time component, = (no, 0, 0, 0). 
Then, upon introducing the convenient notation ^ = (no)^, one can rewrite ([2]) as 

C = ~\e^uF^'' + -L^e^'^'E^.doEki , (3) 

and in particular one can exploit the simplifications provided by spatial isotropy. 

We shall also focus here on this case Ua = (no, 0,0,0), although, unlike previous authors, 
we shall be rather careful (particularly in Subsection 5.1) in assessing the limitations that this 
simplifying assumption introduces. 

From ([3]) one easily obtains modified Maxwell equations for the electric and magnetic fields: 

= djE^ 

= - doE^ - d.e'^'Bi - ^(9o)'S^ (4) 

and 

= dke'^''^^Ei - doB^ 

= djB^. (5) 

Note that this second set of equations is undeformed with respect to the classical (^ = 0) case, 
since it follows only from the antisymmetry of the electromagnetic tensor F^j,. The resulting 
equation of motion for the electric field is: 

= -d^doE' - d\diE^ - d'^Ei) - ^e^'^^idofd^E^. (6) 

For the case of plane waves, Em{x) = Efn{k)e~^''''^p , this equation of motion takes the form: 

= klE^ + k'ihE^ - yEi) - '^(ko)\kme^^^E^) . (7) 

Since the Lagrangian density of Eq. ([3]) is still symmetric under space-rotations there is 
clearly no loss of generality in fixing the direction of propagation of the plane wave in, say, the 
direction z, with wave vector kp = {uj,0,0,p). For such a plane wave propagating along the z 
axis the three equations of motion given by Eq. ([7]) take the form 

= io^E^ 
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= - u^Ei+p^Ei + ^cj^pEa 

Ep 

2i^ 

= - uj^E2+p^E2 - -^uj^pEi. (8) 

Ep 

The first of tliese equations simply states tliat tlie longitudinal component of the electric field 
is still absent (even for ^ 7^ 0), and therefore for a plane wave propagating along the z axis one 
has E3 = 0. The remaining two equations can be compactly reorganized as follows: 

= (-cj2 + {El ± iE2) ± ^to^p (^1 ± iE2) , (9) 

meaning that the right- and left- circularly polarized fields E± = E^e^ ± iE^iy satisfy the 
equation of motion: 

-cj2 +p^± ^to^p] E± = 0. (10) 
J 

So the electric field is subject to birefringence: the left- and right-circular polarized components 
have different dispersion relations: 

cu± ^ p (^1 ± J-p^ , (11) 

where we included only the leading correction in ^/ Ep. The components propagate with different 
group velocities Vg± = ^ 1 ± 2-|-p. 

Assuming that the electric field is originally linearly polarized, then, due to the different 
group velocities of its circularly polarized components, if the time of propagation of the radiation 
is sufficiently long, one might have an observably large difference in the time of arrival of the 
two circularly-polarized components of the field, and the linear polarization would therefore be 
lost [22]. If instead the time of propagation is not sufficient to produce a detectable difference 
in the time of arrival of the two modes with different group velocities, the direction of the 
original linear polarization is rotated (see [27] and references therein). In particular, if the 
wave propagates for a time T, the amount of rotation is: 

a = {lo^ - u+)T = 2^p'T . (12) 
3. Effects on CMB polarization and its povi^er spectrum 

CMB radiation is known to have been generated with some degree of linear polarization (see 
e.g. [28] and references therein), and some linear polarization is observed, so our analysis 
must be performed in the regime governed by Eq. (fT2|) . In light of Eq. (fT2|) . measurements of 
the amount of rotation of the CMB polarization can be interpreted as measurements of the 
parameter ^. To be more precise, one has to consider the fact that the CMB photons propagate 
into an expanding universe, so the energy of the photons is redshifted during the propagation. 
The photon energy redshift dependence is: 

ij{z)=u;o{l + z), (13) 
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where luq is the energy of the photon when it arrives on Earth. So the amount of rotation of 
the polarization plane is, for a CMB photon propagating for a time t from the last scattering 
surface (t = 0): 

a{t)= f2^p{t'fdt\ (14) 
Jo 

The amount of rotation of the polarization of photons propagating from the last scattering 
surface toward us (t = T) can be written in terms of the redshift z as (on the last scattering 
surface z = Z, now z = 0): 

a = 2^pl [\l + z)H-'dz. (15) 
J^P Jo 



Since to a first approximation H = i/o\/^m(l + z)^ + ^a, where Hq ^ 10 ^^s ^ is the Hubble 
parameter (in terms of the reduced Hubble constant Hq = h ■ 100{Km/ s) /Mpc): 

- g r . - (16) 

The formal description of the rotation with angular velocity = ^ = of the CMB 
polarization plane due to any physical effect requires a modification of the Boltzmann equation 
which governs the time evolution of a single Fourier mode of the perturbation of the Stokes 
parameters Q and U. For scalar perturbations the Boltzmann equations become [29| [30]: 



Aq + ikfiAQ 



+ 2ujAi 



Au + ik^Au = - fAu - 2ljAq (17) 

where the dots indicate derivatives with respect to the conformal time t], fi is the cosine of the 
angle between the photon direction and the Fourier wave vector, f = XeneCTd/ao, with Xe the 
ionization fraction, the electron number density, the Thomson cross-section and a the 
scale factor (ao the scale factor today). Sp = At2 + Aq2 — Ago is the source function, A^i and 
Agi indicate, respectively, the i-th moment in the expansion of the temperature perturbation 
At and of the polarization perturbation Ag in Legendre polynomials; P2(At) is the second 
Legendre polynomial. Note that the equation for Au has no generating term, since the U 
component of polarization is not generated by scalar perturbations, but by tensors. 

These coupled differential equations can be separated taking the sum and the difference of 
the first equation with i times the second. Then a formal integration leads to the solution [29] : 

{AQ±tAu){vo) = rfr] f(r/)e''='^(''^''°)-"('')5p(r7)e^'*^''o^'''"^''') (18) 



where r^rj) = J*"' drj'fij]') and rjQ is the conformal time today. 

Since in our case both uj and a are expected to be very small and the visibility function 
f{r])e~'^^^^ peaks tightly at the time of decoupling r] = rjdec, a very good approximation to the 
above solution is given by: 

2 / dr]'uj{r]') j = AQ{r]o) cos (2a) 

2 J dr]'uj{r]')j =-AQ{r]o)sm{2a) (19) 
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where Aq(?7o) = Aq(?7o)|i^=o and a is given by eq. f|T6|) . Taking into account the tensorial 
modes of the perturbations, which generate a primordial U mode, it can be shown that the 
effects of a rotation of the polarization direction are: 

Aq(?7o) = Aq(?7o) cos (2a) + Auivo) sin (2a) 

Au {m)= - Aq (r/o) sin (2a) + Au (vo) cos (2a) . (20) 

This leads to the polarization power spectrum Cf^ ^ J dk [/c^Ax(^7o)Ay(r7o)], X,Y = 
T,E,B: 

Cf^ = Cf^ cos^ (2a) + Cr sin' (2a) 
^ = Cf^ sin' (2a) + ^ cos' (2a) 

Cr = I {Cr-Cr) sin {4a) 

Cf^ = Cj^ cos (2a) 

Cj^ =Q^£^sin(2a) (21) 

We checked numerically that evaluating the power spectra given by equation ( 1211) instead of 
the exact formula in equation (fT8|) (and the corresponding one for tensorial modes), leads to 
an error on each Cf^ that, for J^^"" drj'uj{rj') up to six degrees, is much less than 1%. 

4. Data Analysis 

In this section we discuss the current constraints on a and on ^ from the most recent CMB 
polarization datasets coming from the WMAP five-year mission [HI [31] and the BOOMERanG 
polarization flight in 2003 [32]. We follow the standard approach in the literature by making 
use of the the publicly available Markov Chain Monte Carlo package cosmomc [33] with a 
convergence diagnostics done through the Gelman and Rubin statistic. 

In our analysis we consider the complete set of power spectra (12T|) . taking into account 
also TB and EB cross-correlation power spectra, which are usually set to zero in the standard 
CMB calculations. We sample the following eight-dimensional set of cosmological parameters, 
adopting fiat priors on them: the physical baryon and Cold Dark Matter densities, ujb = VL^h? 
and ujc = fich'^, the ratio of the sound horizon to the angular diameter distance at decoupling, 
Og, the scalar spectral index Ug, the overall normalization of the spectrum Ag at k = 0.05 Mpc~^, 
the optical depth to reionization, r, the tensor to scalar ratio of the primordial spectrum, r 
and, finally, the rotation of the power spectrum discussed above, a. We derive the value of ^ for 
each Markov chain step using f|T6|) . so including the dependence on the geometrical parameters 
of the universe. Simultaneously we also use a top-hat prior on the age of the Universe as 
lOGyr < to < 20 Gyr. Furthermore, we consider purely adiabatic initial conditions, we impose 
flatness, we treat the dark energy component as a cosmological constant, and we include the 
weak lensing effect in the CMB spectra computation. 

In Table [T] we show the current constraints on a and ^ from current experiments. WMAP 
5- year data provide the constraint of a = (—1.6 ± 2.1)° (consistent with previous analyses 
[T8l [H]), i.e. with no indication for a different from zero, while BOOMERanG data give 
a = (—5.2 ±4.0)°, with a ~ 1.3cr hint for a < 0. Assuming a energy-independent (a = ao), the 
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joint estimation from BOOMERanG and WMAP5 data gives ao = (— 2.4± 1.9)°, which is also 
consistent with [2]. As regards to the ^ parameter, from WMAP 5-year data we constrain it 
to be ^ = -0.09 ±0.12, while with BOOMERanG datcHwe obtain ^ = -0.123 ±0.096. Since ^ 
is energy- independent, we can also make a combined analysis of the two datasets, which results 
in the estimation of ^ = —0.110 ± 0.075. 



Experiment 


a ± cr{a) 


e±^(0 


WMAP (94 GHz) 


-1.6 ± 2.1 


-0.09 ± 0.12 


BOOMERanG (145 GHz) 


-5.2 ± 4.0 


-0.123 ± 0.096 


WMAP±BOOMERanG 




-0.110 ± 0.075 



Table 1. Mean values and lu error on a(in degrees) and ^ for WMAP, BOOMERanG and 
WMAP+BOOMERanG 

We also analyze the ability of future experiments in constraining a and ^. We created 
several mock datasets with noise properties consistent with the PLANCK satellite [3lj, the 
Spider balloon borne experiment [35] and the EPIC satellite [36], assuming, as fiducial model, 
the best fit of WMAP with a = (see TableE] for the experimental specifications used in the 
analysis) . 

We analyze these datasets with the full sky exact likelihood routine [37j in the cosmomc 
code including the /f^^ factor to reduce the degrees of freedom, ignoring correlations between 
different multipoles i. Firstly, we consider the specifications for the PLANCK satellite for the 
channels at 70, 100, 143 and 217 GHz respectively. Clearly the presence of different channels 
in the PLANCK experiment has the main goal of foreground removal. However we report the 
sensitivity on a and ^ for each channel considering the possibility of a frequency dependence of 
the cosmic signal. As we can see in Table [3] the PLANCK channel with the highest sensitivity 
on a is the 143 GHz channel, able to constrain a at the level of 0.07° and ^ at level of 0.0017. 
Lower sensitivities on a will be reached by other "high frequency" channels but, because of 
the energy dependence of a, the 217 GHz channel is the most sensitive to ^, since it allows 
to constrain this parameter at the level of 0.0010. The LFI 70 GHz channel can provide a 
constraint of the order of 0.64° for a and 0.06 for ^. It is interesting to observe that combining 
the three HFI channels (100, 143, 217 GHz) it is possible to improve the sensitivity on the ^ 
parameter, up to 8.5 ■ 10~^. 

Before continuing it is important to bear in mind that any systematic present in the 
calibration of the polarization angle will be reflected in a spurious a and, as a consequence, a 
spurious ^. The current HFI and LFI channels are expected to be calibrated at a level of about 
~ 1° thanks to measurements of polarized emission from the Crab nebula. The expected values 
reported in the Table are therefore not taking in to account this possible systematic. However, 
a multi-frequency analysis could permit to distinguish a systematic effect due to miscalibration 
from a true birefringence effect, since the first is energy independent, while the second has a 
specific energy dependence. 

^ For the analysis using the BOOMERanG dataset, we fixed the optical depth r to the fiducial value of r = 0.09. 
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Experiment 


Channel 


FWHM 


AT/T 


PLANCK 


70 


14' 


4.7 


ftkit — 0.85 


100 


10' 


2.5 




143 


7.1' 


2.2 




217 


5' 


4.8 


Spider 


145 


40' 


0.30 


f J. — 0.50 










70 
/ u 


1 9' 


n nfifi 

U.UUL) 




100 


8.4' 


0.066 




150 


5.6' 


0.084 




220 


3.8' 


0.17 



fsky — 0.85 



Table 2. PLANCK, Spider and EPIC experimental specifications. Channel frequency is given 
in GHz, FWHM in arcmimitcs and noise per pixel in 10^^ for the Stokes I parameter; the 
corresponding sensitivities for the Stokes Q and U parameters are related to this by a factor of 
^/2. 



Exi)eriniont 


Cliannel 


a{a) 




(T 


(0 




PLANCK 


70 


0.64 




6.0- 


10- 


-2 




100 


0.14 




6.5- 


10- 


-3 




143 


0.073 




1.7- 


10- 


-3 




217 


0.10 




1.0- 


10- 


-3 




100+143+217 






8.5- 


10- 


-4 


Spider 


145 


0.27 




6.1 - 


10- 


-3 


EPIC 


70 


2.1 • 10" 


-3 


1.9- 


10- 


-4 




100 


1.8-10- 


-3 


7.8- 


10- 


-5 




150 


1.5-10- 


-3 


2.9- 


10- 


-5 




220 


1.2 - 10- 


-3 


1.1 - 


10- 


-5 




70+100+150+220 






1.0- 


10- 


-5 


CVL 


150 


6.1 ■ 10- 


-4 


1.3- 


10- 


-5 




217 


6.1 - 10- 


-4 


6.1 - 


10- 


-6 



Table 3. Expected Ict error for PLANCK 70, 100, 143, 217 GHz, Spider 145 GHz, EPIC 
70, 100, 150, 220 GHz and two ideal CVL experiment at 150 GHz and 217 GHz on a (in degrees) 
and ^. 

It is interesting to extend the forecast to other future experiments as Spider and EPIC. 
Spider, a balloon borne experiment that will fly over Antarctica in 2012, will constrain variation 
of a at level of 0.27" and of ^ at level of 6.1 ■ 10^^, with a sensitivity competitive with PLANCK. 
The EPIC satellite proposal, or an equivalent next generation CMBpol satellite mission, can 
detect deviations as small as 0.0012° for a and 1.0-10^^ for S,-, providing a dramatic improvement 
respect to PLANCK and Spider. Again, the channel at higher frequency provides the best 
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constrain. Finally we consider an ideal experiment, cosmic variance limited in anisotropy 
and polarization measurements. This experiment provides a sort of fundamental limit to the 
precision that can be achieved. We found that in this, ideal, experiment, angles as small as 
a = 0.0006° could be, in principle, measured, and we observe in Table Othat this would provide 
sensitivity to = 1.3 ■ 10~^, for an experiment at 150, or even ^ = 6.1 ■ 10~^, for an experiment 
at 217 GHz. 

5. Comparison with other results 

5.1. Comparison with limits previously obtained using observations in astrophysics 

While ours is the first explicit work aimed at comparing the predictions of the Myers-Pospelov 
framework to CMB data, there have been a few other attempts, such as the ones reported in 
Refs. [22l [71 [231 [25] , to place limits on the parameter ^ using certain observations in astrophysics. 
The outcome of some of these analyses was simply described in terms of limits on the parameter 
^, as if they were absolute limits, but we must here stress that, in light of the very explicit 
frame dependence of the Myers-Pospelov model, one of course cannot establish absolute limits 
on parameters. From the analysis of experimental data collected in a certain frame one can 
definitely obtain even very robust limits on the values of some of the components of the four- 
vector Ua (see Eq. ([2])) in that frame, but it is highly nontrivial to then qualify such results 
from a frame-independent perspective. 

All these previous studies have, like the one we are here reporting, focused on the 
implications of a time component for the Myers-Pospelov symmetry-breaking four- vector n^, 
a restriction which is introduced very explicitly by assuming = (no? 0,0,0). But of course 
the quantum-gravity arguments [21 [3l [2T] that motivated the study of the Myers-Pospelov 
framework do not provide any indication that the four-vector should take this form in any 
specific frame {ua may well not even be time-like), and even if it took this form in a certain 
class of frames it would then definitely have different form in other frames, transforming like a 
four- vector from frame to frame. 

Constraints on all four components of Ua in one frame could be converted, through 
appropriate Lorentz transformations, into analogous (but possibly very different in magnitude) 
constraints applicable in another frame. But bounds established exclusively for the time 
component of are of very limited applicability in other frames. 

This is particularly significant for our results since some analyses of data in astrophysics 
can lead to very stringent bounds on ^. In one recent such study [25] it was even argued that 
observations of polarized radiation from the Crab Nebula can be used to obtain the impressive 
bound 1^1 < 10^^, which would amount to the constraint \no\ < lO^"' on the time component of 
the four- vector n^. Clearly these results obtained in astrophysics start to provide us an intuition 
that large values of uq are disfavored, but at the same time it should be noticed that it is not 
easy to convert such results into intelligible constraints applicable in reference frames that are 
different from the "laboratory frame" where the bound was derived. In principle one could even 
contemplate the possibility that in such a laboratory frame one has, say, Ua = (0, 10^, 10^, 10^), 
which would be compatible with any upper bound on no and would still be meaningful from a 
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quantum-gravity perspective (since, as mentioned in Section 2, components of with values of 
order 10^ could be expected in rather popular scenarios estimating the scale of quantum-gravity 
as roughly given by the grand-unification scale). And if indeed Ua = (0, 10^, 10^, 10^) in some 
laboratory frame then a boost of unimpressive magnitude, with, say, /? ~ 10~^, could take us 
from that laboratory frame to a frame where the time component of is actually of order 1. 

In light of the rather significant frame dependence implied by these considerations it 
is rather clear that for this research program it is a top priority to move on to bounding 
experimentally all four components of Ua- And for works that focus on the time component it 
appears that it is rather advantageous to focus on studies of the CMB since one can at least 
rather easily combine different sets of data in a meaningful way, by describing them all in the 
natural reference frame of the CMB. 

Also significant from the perspective of the frame dependence of the Myers-Pospelov 
framework is the fact that the results of our analysis, while being inspired by this framework, 
depend most strongly on the assumption of a characteristic (quadratic) energy dependence of a, 
and not very sensitively on the details of the Myers-Pospelov framework with = (no, 0, 0, 0). 
Such an energy dependence is a generic feature of any quantum-gravity-inspired description of 
birefringence, since it is obtained using only dimensional analysis and the fact that the effects 
must disappear in the Ep ^ oo limit. One should therefore expect a similar energy dependence 
also in the Myers-Pospelov framework with generic and in models where the breakdown of 
Lorentz symmetry is governed by, say, a two-index tensor. 



5.2. Comparison with previous CMB limits on energy- dependent birefringence 

While ours is the first study to discuss a robust link from the parameters of a model that has 
been extensively studied in the quantum-gravity literature and CMB polarization data, from 
a broader phenomenological perspective the possibility of a quadratic energy dependence of a, 
which is indeed the main characteristic of the effects we studied, had already been considered 
in Ref. [20]. In order to compare of our results to the ones of Ref. [20] we must establish a 
link between the parameter ^ of the (isotropic) Myers-Pospelov model and the multi-parameter 
description of birefringence effects for CMB polarization data analyses adopted in Ref. [20] . This 
can be done by describing the rotation of linear polarization of CMB radiation codified in our 
Eq. (fT6|) in terms of rotations of the Stokes vector s = (Q, U, V)'^ , since it is at the level of these 
rotations of the Stokes vector that Ref. [20] introduced its multi-parameter phenomenological 
picture. 

Let us start by noticing that the variation of the Stokes vector due to birefringence (not 
necessarily between circularly polarized components) can be written as |20j : 

ds 

— = 2cuC X s (22) 
dt 

where ( is proportional to the Stokes vector of the faster of the two modes subject to 
birefringence. 

The Stokes vector ( can be written in the spin- weighted basis, (sw = (C(+2), C(o)5 C(-2)) = 
(C^ — zC^, + iC^)^) ^-iid each component of (sw can be expanded in spin-weighted spherical 
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harmonics: 

C(±2) = {k{E)lm ± ik(^B)lm) ±2^Zm, C(0) = k(^YyrnoYlm- (23) 



Im Im 

Each coefficient k(y)im and, respectively, k(^E,B)im in Eq. (123|) can be expressed as a 
ombination of only odd and, respectively, even powers of uj. The parameters k^^^i^ and 



^(E,B)im ^s^*^ EH] are the coefficients of the expansion: = Y^dodd^ ^{v)im^ k(^E,B)im = 

V ,.d-4l.id) 
Z^deven^ '^{E,B)lm- 

In the case of our interest, in which the rotation of the Stokes parameter is due to a 
birefringence between circularly polarized components of the radiation, taking conventionally 
as the fastest mode the right-circularly polarized one, one finds that the vector C, is proportional 
to (0,1,0) in the spin-weighted basis. The variation of the Stokes vector is, in the (Q, Vy 
basis: 

^ — Qo sin(a) -|- f/o cos(a) ^ 



ds 
di 



—Qo cos(q;) — Uq sin(Q;) 
V y 

where a{t) is given by Eq. f|T4|) : 

Writing this equation in the spin-weighted basis, evaluating it at time t = T (today) and 
comparing with Eq. fl22l) . one finds C{o) = ~e~Po- ^'^ coefficients of the expansion 

of C(±2) are null, while for ^(o) the only nonzero coefficient is k(v)oo = ~e^PoV^, because 
C(o) doesn't depend on the direction of observation in the sky. Since k^v)oo is linear in the 
energy of the photon, there is only one nonzero coefficient in the expansion of k(^v)oo, and this 
coefficient, which was denoted by fc[y)oo in Ref. [20], is therefore the only parameter in the 
parameterization introduced in Ref. pO] that would reflect the presence of a (time-component- 
only) Myers-Pospelov term in the Lagrangian density of electrodynamics. The comparison 
between our results and the ones of Ref. [2(^ should therefore be based on the relationship 
between our parameter ^ and the parameter fc(\/-)oo' which is 

For this parameter Ref. j20j arrives at an estimate fc[y)oo — (3 ± 2) x lO^^^GeV^^^, using 
exclusively BOOMERanG data and without performing a Markov chain analysis, but rather 
fixing the cosmological model to the WMAP best ffi model and neglecting therefore possible 
correlations between the cosmological parameters El- In light of ( l25i) our result ^ ^ —0.110 ± 
0.076 amounts to k^fy^^^ = (3.2 ±2.1) x 10~'^^GeV~^ . This is fully consistent with the findings 
of Ref. [20], and we would like to argue that our result should now be viewed as the benchmark 



+ An effect of this could be the fact that another parameter studied by [50] , fc|y>)Qo , is estimated also from the 
BOOMERanG dataset to be fc|y\oo ^ (^'^ ± 7) x lO^^'^GeV , which corresponds to a rotation angle a = 12 ± 7 



^(V)00 

degrees, while the estimate of ^[y^gg = (3 ± 2) x lO^^'^GeV^^ from the same dataset corresponds to a rotation 
angle a = 4 ± 3 degrees. 
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for the limits on /C(y)Qo, both because we rehed on a more sizeable data sample, combining 
BOOMERanG and WMAP data, and because of the marginalization over the remaining cosmic 
parameters. 

6. Conclusions and Outlook 

The main objective of our analysis was to establish robustly that cosmology and particularly 
CMB data can provide constraints on effects that are of interest from a quantum-gravity 
perspective and are introduced at the Planck scale. This possibility has been contemplated 
in several previous studies but never in a way that would transparently expose sensitivity to 
effects introduced truly at the Planck scale. We feel that we have fully achieved our goal by 
relying on the Myers- Pospelov model, and on the rather large quantum-gravity literature on this 
model, which also provided a clear target for the description of effects introduced at the Planck 
scale. We therefore hope that our analysis will motivate a more intense program of studies of 
quantum-gravity effects in cosmology, not necessarily focused on in- vacuo birefringence and/or 
the Myers-Pospelov model. 

It is intriguing that our estimate of the Myers-Pospelov ^ parameter, ^ = —0.110 ± 0.076, 
provides a faint indication for a nonzero effect. In this respect it is certainly valuable that, as 
discussed in Section 4, future experiments will be able to measure CMB radiation polarization 
with higher sensitivities, and we find that this should lead to an improvement of three orders 
of magnitude in the constraints obtained following our strategy of analysis ~ 10~^). 

As stressed in Subsection 5.1, in light of the rather significant frame dependence introduced 
by the Myers-Pospelov setup, this type of sensitivities obtainable from CMB-data analyses are 
still rather meaningful, in spite of the apparently ultra stringent constraints that have been 
obtained on the parameter ^ using astrophysical observations. We also showed that CMB data 
can be used to determine the sign of ^, a possibility which was never exposed in previous works 
relying on astrophysics and that may prove challenging in those contexts. 

From the methodological aspect, while in the past the birefringence of CMB photons has 
been widely studied, and several bounds have already been provided on the birefringence angle 
a, here wc pointed that in the literature the possible energy dependence of the effect has not 
been adequately contemplated, leading to possible misleading conclusions when averaging on 
different photon frequencies. We showed how exploiting the availability of data from different 
channels of frequency of the photons improves the sensibility on the birefringence parameter 
^ and leads to a slightly stronger indication for a nonzero effect, than the results obtained 
estimating the birefringence angle a with an average on frequencies. This should be taken 
into account in planning future CMB experiments, and may turn out to be another reason for 
combining cosmological and astrophysical investigations of the Myers-Pospelov framework. 

Here, we have analyzed a uniform, isotropic rotation to the EM polarization vectors. But, 
as stressed in Subsection 5.1, the quantum-gravity motivation for investigating the Myers- 
Pospelov framework and the intrinsic structure of the Planck-scale correction term added to the 
Lagrangian density of electrodynamics provides a strong invitation to contemplate also effects 
resulting in anisotropic rotation of polarization. In an upcoming paper we plan to return to 
this issue and address the extent to which CMB polarization measurements can constrain such 
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an anisotropic model. 
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